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Abstract. We detect seven distinct absorbing systems in the high-resolution X-ray spectrum of the Seyfert 1 galaxy IC 4329A, 
taken with XMM-Newton. Firstly we detect absorption due to cold gas in our own Galaxy and warm gas in the Galactic halo or 
the Local Group. This local warm gas is only detected through O VII absorption, from which we deduce a temperature between 
0.03 and 0.2 keV. In IC 4329A we detect absorption from the host galaxy as well as from a warm absorber, close to the nucleus, 
which has 4 components. The absorption from the host galaxy is well modeled by neutral material. The warm absorber detected 
in IC 4329A is photoionized and has an ionization range between log f = — 1.37 and log £ = 2.7. A broad excess is measured 
at the O VIII Lya and N VII Lya emission lines, which can be modeled by either disklines or multiple Gaussians. From the 
lightcurve we find that the source changed luminosity by about 20 % over the 140 ks observation, while the spectral shape, i.e. 
the softness ratio did not vary. In the EPIC spectra a narrow Fe Ka and Fe XXVI Lya emission line are detected. The narrowness 
of the Fe Ka line and the fact that there is no evidence for flux variability between different observations leads us to conclude 
that the Fe Ka line is formed at a large distance from the central black hole. 
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•• • 1. Introduction in NGC 3227 and Ark 564. These similarities are the high in- 

, ' clination angle and the reddening determined between the far 

k> IC 4329A is one of the X-ray brightest Seyfert 1 galaxies in the uv and optical band Crenshaw & Kraemer J200T1 predict that 

r £ ■ sky. The spiral galaxy is seen almost edge on and is the closest only lowly ionized absorption shou i d be detected, namely the 

C3 ; companion to the elliptical galaxy IC 4329, the central galaxy dust that causefj the reddening observed in the optical and the 

in a cluster. IC 4329A has rather extreme properties, such as the jj y band 
full width at zero intensity (FW0I) of the Ha line which equals 

13,000 km s- 1 (Disney, 1973}. On the basis of the reddening T !^ shift ° f IC 4329A is 01605 ± 00005 ( Willmer 
observed in the optical, A v = 2.5 - 4.8 magnitudes, Wilson & et aL M)- However, Wilson & Penston {TW^ note that the 
Pension (fl979l have classified IC 4329A as the nearest quasar, °P tlcal lines are double P eaked ' Wlth the second P eak at a mean 
with an absolute magnitude between -23.0 and -25.3. In the redshlft of 0224 ± 0001 ■ Throughout our analysis we fixed 
optical spectra prominent Na I D, Ca II and weaker Mg I ab- the redshlft to 2 = 0-01605, as we did not detect any evidence 
sorption lines are observed, indicating the presence of dust, for a component at z = 0.0224. All spectra were corrected for 
consistent with the dust lane observed in the equatorial plane a Galactic hydrogen absorption column density N H = 4.55 x 
of the galaxy (Wilson & Penston[l973- The optical spectrum 10 m 2 ( Elvis ' Lockman & Wilkes[l98l. 
shows that the absorber is neutral or very lowly ionized. The Earlier X-ray spectra contained a reflection component 
observed line emission is due to lowly ionized material with (Miyoshi et al. U988t and a moderately broadened Fe Ka emis- 
strong [O III] (Wilson & Penston fT979l . Crenshaw & Kraemer sion line (Done et al. l2000> . In an earlier 12 ks XMM-Newton 
(2001 1 classified IC 4329A as a dusty luke-warm absorber on observation Gondoin et al. J2001t detected only a narrow Fe 
the basis of similarities with the known absorbers of this type Ka line, with a width a = 0.01 ± 0.05 keV. The earlier Chandra 
HEG spectrum shows a complex double line at the Fe Ka en- 
Send offprint requests to: K.C. Steenbrugge ergy, which can be fit either by two Gaussians, two disklines, or 
Correspondence to: K.C.Steenbrugge@sron.nl a single diskline (McKernan & Yaqoob 2004 1. Gondoin et al. 
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Table 1. Instrumental set-up and exposure time used for the 
spectral analysis. 



Instrument 


exposure (s) 


observational mode 


RGS 1 and 2 


136017 


spectroscopy 


pn 


44363 


Small window with thin filter 


MOS 1 


97857 


PrimePartialW3 with medium filter 


MOS2 


94709 


Small window with thin filter 



(2001 1 fitted the soft X-ray part of the spectrum with a set of 
absorption edges of O I, O VI, O VIII and N VII. They noted the 
lack of an O VII absorption edge and concluded that the warm 
absorber must have at least 2 phases. 

In this paper we analyze a deep 140 ks XMM-Newton ob- 
servation of IC 4329A. In Sect. 2 we describe the observation 
and data reduction. The variability is analyzed in Sect. 3. In 
Sect. 4 the spectral analysis is detailed, while Sect. 5 details 
the analysis of the absorbers present in the spectrum. In Sect. 

6 we discuss our results and the conclusions are summarized 
in Sect. 7. Power density spectra analysis as well as possible 
variability in the Fe Ka line will be described by Markowitz 
et al. (20061 once the complete RXTE long term monitoring 
campaign is completed. 

2. Observation and data reduction 

IC 4329A was observed on August 6 and 7, 2003, for a to- 
tal exposure time of 136 ks. Throughout the observation there 
were periods of high background, which were removed for the 
EPIC data using the Good Time Interval (GTI) produced with 
the online Processing Pipeline Subsystem PPS data files. For 
the RGS spectra we use the whole exposure time, to maxi- 
mize the statistical quality of the data. For the data reduction 
we used the PPS data products, downloaded from the XMM- 
Newton website. The data were reduced with the SAS software 
version 5.4.1. Tabled lists the instrumental set-up and the net 
exposure time used in the further spectral analysis. For the time 
variability analysis we used the full 140 ks pn observation. The 
pn, MOS 1 and MOS 2 spectra were fitted between 0.4 - 12 
keV, the first order RGS spectra were originally fitted between 

7 — 38 A, the second order spectra between 7 — 19 A. Due 
to the low count rate, we limited our analysis to the first order 
spectra of RGS, and the 7 to 26 A wavelength range. RGS 1 
and 2 were averaged in the presented plots for clarity, but fit 
as separate spectra. For both RGS instruments the count rate 
ranged between 30 and 250 counts per bin for the part of the 
spectrum up to 26 A, allowing the use of x 2 statistics. The RGS 
and MOS spectra were binned by a factor of 3, for the RGS this 
corresponds to about half the FWHM resolution of 0.07 A or 
a bin size of 0.03 A. The MOS cameras have a FWHM reso- 
lution of 50 eV at 0.4 keV and 180 eV at 10 keV. The pn data 
were binned by a factor of 5, and have a FWHM resolution 
ranging between 90 eV and 185 eV at 1 keV and 10 keV re- 
spectively. All the spectra were analyzed using the SPEX soft- 
ware (Kaastra et al. 2002b I. All quoted errors are rms errors, 
i.e. Ax 2 = 2. 



absorbed spectrum on IC 4329A 




Time since start (in ks) 



Fig. 1. The background subtracted lightcurve for the pn instru- 
ment, using the complete 140 ks observation. 
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Fig. 2. Lightcurves for the pn instrument. Top: the hard band 
between 2—10 keV; middle: the soft band between 0.4 — 1 
keV. Bottom: the softness ratio: the 0.4 — 1 keV count rate 
divided by the 2—10 keV count rate. All the data were binned 
in 1000 s bins. 



3. Time variability 

We constructed the light curve using the pn in the 0.4 — 10 
keV band (see Fig.0. Time is measured from the start of the 
observation, at the 6 th of August 2003 at 6 h 13 m 47 s UT. Fig.|2] 
shows the lightcurves in the hard (2—10 keV) band, the soft 
(0.4 — 1 keV) band and the softness ratio determined from 
these two bands. 

The count rate of the pn varied between 20 counts s _1 at 
the start of the observation and 16 counts s _1 at the end of the 
observation. The most noticeable part of the light curve is the 
steady decline lasting 25 ks near the end of the observation, 
after which it levels off for the last 3 ks. Throughout the obser- 
vation several gradual luminosity changes are observed. 
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Fig. 3. RGS spectrum of IC 4329A with the strongest absorp- 
tion edges identified. RGS 1 and RGS 2 are averaged for clar- 
ity. Note that the edges are convolved with the resolution of the 
RGS and are heavily blended, thus they are not sharp. 



We identify in Fig. several parts in the light curve where 
the total luminosity decreases or increases nearly monotoni- 
cally. The smallest duration for a 5 % luminosity variation is 
~ 10 ks. The characteristic timescale, r = |dt/dlni |, for these 
luminosity variations is between 100 and 300 ks. The r.m.s. 
variance measured as the intrinsic spread for this dataset is 5 %. 
From the softness ratio (see Fig.|2ji no time lag can be discerned 
and there is no spectral change with increasing luminosity ob- 
served. 



4. Spectral analysis 

4.1. Continuum 

The spectrum is heavily absorbed, with several deep edge-like 
structures (see Fig.|3]for RGS and|4]for EPIC-pn). Below en- 
ergies of about 1.3 keV (i.e. above ~ 10 A) we need a com- 
plex absorption model (see Sect. |5j, a power law continuum 
and a soft excess which we successfully model with a high 
temperature modified black body. This is a black body com- 
ponent modified by coherent Compton scattering, for more de- 



tail see Kaastra & Barr (1989 1. Above 1.3 keV the spectrum 
is well described by a power-law plus two Gaussians for Fe 
Ka and Fe XXVI Lya emission line probably blended with the 
Fe K/3 line (see Sect. 14. 2> . corrected for Galactic absorption 
(see Fig.|4]and Fig[6}. We modeled the Galactic absorption to- 
ward IC 4329A with an absorber in collisional ionization equi- 
librium, namely the hot model in SPEX. The hot model de- 
termines the transmission for a certain temperature and total 
hydrogen column density. This model includes more accurate 
atomic data for the oxygen edge and includes the neutral oxy- 
gen absorption lines. We froze the temperature to 0.5 eV (re- 
sulting in an almost neutral gas) and the column density to JVh 
= 4.55 x 10 24 rar 2 (Elvis et al. IT9891 We froze the abun- 
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Fig. 4. The pn spectrum of IC 4329A. Note the strong absorp- 
tion at lower energies (see Fig.[6]for the residuals to a power- 
law fit) and the strong Fe Ka emission line. 

Table 2. The best fit continuum parameters, for the fit including 
the absorbers detailed in Sect. |5] fitting the pn and RGS 1 and 
RGS 2 simultaneously. 



pl 


norm." 


2.90 x 10 52 ph s- 1 keV' 1 




lum. a ' b 


1.17 x 10 37 W 






1.709 


mbb 


norm." 


7.2 x 10 31 m 1/2 




kT a 


473 eV 



The systematic errors dominate over the statistical errors, which are 
on the order of 1%. 
b in the 2—10 keV band, calculated from the measured 
normalization, photon index. 



In our spectral analysis of the low energy part of the spec- 
trum we focus on the pn and RGS spectra. The pn spectrum 
is used to determine the continuum and broad emission line 
parameters, while the RGS spectra are used for fitting the dif- 
ferent absorbers in detail. We did not fit the MOS 1 and MOS 
2 spectra simultaneously with the pn and RGS spectra as the ~ 
10% systematic uncertainties in the effective areas between the 
different instruments influence the fit noticeably. The best fit 
absorption model obtained from the pn and RGS spectra alone 
matches the MOS spectra rather well, if we allow for different 
continuum parameters. Table[2]lists the best continuum param- 
eters for the simultaneous pn, RGS 1 and RGS 2 fit. We did 
not include a continuum reflection component as detected by 



dances to solar and the velocity broadening, the Gaussian a to 
the standard 100 km s _1 . 



Gondoin et al. (2001 1 from the BeppoSAX data in our spectral 
analysis. The XMM-Newton spectra only extend to about 12 
keV, where the continuum reflection component is still negligi- 
ble (see Fig. [6). 

4.2. Fe Ka and Fe xxvi Lya line 

A narrow Fe Ka line is clearly detected in all three EPIC spec- 
tra. There is a second, weaker line detected at about 6.9 keV, 



4 



K. C. Steenbrugge et al.: Heavily absorbed spectrum on IC 4329A 



Table 3. The best fit parameters for the Fe Ka and Fe XXVI 
Lya emission lines for the pn and MOS 1 and MOS 2 spec- 
tra, using Gaussian profiles. As the line at ~6.9 keV is rather 
weak, we fitted MOS 1 and MOS 2 simultaneously. The ener- 
gies quoted are for the restframe of IC 4329A. The errors on 
the line energies are the statistical uncertainties since these are 
larger than the systematic uncertainties in the energy scale. All 
quoted results have been corrected for the redshift of IC 4329A. 





pn 


MOS1 


MOS2 


norm™ 


9.5 ± 0.9 


7.6 ± 0.9 


6.4 ± 0.8 


flux 6 


0.85 ± 0.08 


0.66 ± 0.08 


0.58 ± 0.07 


E (keV) 


6.45 ± 0.01 


6.42 ± 0.01 


6.40 ± 0.01 


FWHM (eV) 


240 ± 50 


110 ±50 


< 100 


EW (eV) 


87 ± 8 


68 ± 8 


64 ±8 


norm™ 


3.0 ± 0.8 


3.2 ± 0.9 


flux 6 


0.27 ± 0.07 


0.28 ± 0.08 


E (keV) 


6.92 ± 0.05 


6.89 ± 0.05 


FWHM (eV) 


300 ± 100 


300 ± 100 


EW (eV) 


30 ±8 


36 ± 10 



™ in x 10 49 phs" 1 . 
6 inph m -2 s _1 



which we identify as a blend of the 6.95 keV and 6.97 keV 
Fe XXVI Lya lines and probably blended with the Fe K/3 line 
at 7.06 keV (Palmeri et al. l2003l . In Fig.|5]the pn, MOS 1 and 
MOS 2 spectra around the Fe Ka line are plotted, in Fig.[6]the 
residuals to a power-law fit of the pn spectrum is shown. 

We tested the accuracy of the energy scale for all instru- 
ments by looking at instrumental background lines, to avoid 
any gain problem. For the MOS cameras the outer CCDs are 
exposed, and we used the whole exposed area to find the in- 
strumental background lines. For MOS 1 the strongest instru- 
mental line, Al at 1.4866 keV is detected at the expected en- 
ergy. For MOS 2 all the background lines, with the exception 
of the weak Ni line at 7.4724 keV are detected at the expected 
energy. Therefore the energy scale of the MOS spectra should 
be correct. For pn we are limited to CCD 4; however, we do not 
detect any of the the instrumental background lines with suffi- 
cient significance to be able to verify the energy scale. There 
can thus still be a small gain problem for the pn. 

Table [3]lists the best fit values fitting Gaussian profiles to 
both iron emission lines in the pn and the MOS spectra. The 
full width half maximum (FWHM) and energy of the Fe Ka 
line are consistent within the instrumental uncertainties. For the 
line at 6.9 keV data for the pn and the MOS instruments are 
consistent. Interestingly, we do not detect the FeXXV triplet 
around 6.70 keV. 

Gondoin et al. (2001 1 find an absorption edge at 7.1 keV 
with an depth of 0.03 ± 0.02. There is a 1 a jump in our data at 
7.1 keV, corresponding to an optical depth of < 0.04. However, 
we do not detect sufficient iron L line absorption between 15.5 
A and 16.5 A and edge absorption between 8 A and 10 A in the 
RGS spectra to explain such an edge. As a result we decided 
against the interpretation that this feature is an edge. McKernan 
& Yaqoob (2004i note two emission features in the Chandra 




5.5 6 6.5 7 7.5 



Energy (keV) 

Fig. 5. Detail of the pn, MOS 1 and MOS 2 spectra showing the 
Fe Ka and Fe XXVI Lya emission lines. The model through the 
data is given in Tabled The MOS 2 data were shifted by adding 
0.006 counts s _1 for clarity. 
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Fig. 6. A power-law plus Galactic absorption fit to the pn spec- 
trum, rebinned by a factor of 9. Note the Fe Ka and Fe XXVI 
Lya emission lines, as well as the absorption at lower energies. 
The feature at 2 keV is an artifact of the absorption spectrum 
and the soft excess. Note that the reflection component has at 
most a 1 a significance. 

HEG spectrum, which they model with either two Gaussians or 
two disklines. In the model with two Gaussians they obtain for 
the Fe Ka line an energy of 6.30 ± 0.08 keV, for the Fe XXVI 
Lya line an energy of 6.91 ± 0.04 keV. Both are consistent 
with our results, as are the measured EWs and intensities. We 
also tested a double diskline model (Laor 1991 1, which gives a 
similarly good fit, for two more free parameters as compared to 
the fit with two Gaussians. For pn we find an inclination angle 
i < 17° and an emissivity slope q < 1.4, where for large radii 
J(r) ~ r~ q . This is consistent with the parameters determined 
by McKernan & Yaqoob (2004i who find an inclination angle 
of 12° ± 1 ° and a flat emissivity slope. For the MOS spectra we 
find similar, but less constrained upper limits. The diskline pro- 
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Fig. 7. Models for the two broadened emission features de- 
tected in the RGS spectrum (see Sect 15. 4> and the two emis- 
sion features detected in the pn spectrum. The thin solid line 
represents the best fit for the relativistically broadened lines 
detected in the RGS spectra, with q = 1.9 and i = 80.9 (see 
Table The dashed line represents the model for the emis- 
sion features adopting the parameters found for the relativisti- 
cally broadened Fe Ka and Fe XXVI emission lines: emissivity 
slope q = 1 and i = 12°, consistent with the values obtained 
by McKernan & Yaqoob (20041. Note that this profile is very 
similar to a Gaussian profile. 

file for i = 12° and q = 1 is shown in Fig.[7](dashed line). Note 
that this profile is very similar to a Gaussian line profile. As a 
result, the pn and MOS spectra do not allow us to distinguish 
between a model with two disklines or two Gaussian emission 
lines. As we cannot constrain the diskline parameters, we use 
the model with two Gaussians in the following data analysis. 

We re-analyzed the pn spectrum for the 2001 XMM- 
Newton observation of IC 4329A. We used S AS version 6.0 and 
the PPS data files for the data reduction. The data are noisy (to- 
tal exposure time of 9.7 ks), leading to large error bars or upper 
limits for the Fe Ka parameters, even if we fit the emission line 
with a Gaussian. In Table0]we compare the different literature 
values for the Fe Ka emission line. Note that the flux of the Fe 
Ka line is consistent with a constant value since August 1997, 
although there is a rather large spread in FWHM. However, the 
error bars on the flux values are large, and thus we can not 
exclude variability. The large spread in FWHM is partially a 
result from differing spectral resolution of the different instru- 
ments. An indication of the uncertainty in absolute calibration 
between the different instruments are the differing parameters 
for the simultaneous ASCA and RXTE observations (Done et 
al. l2000l 

5. Absorbers 

The RGS spectra are dominated by absorption edges (see 
Fig. |3}, as was also concluded by Gondoin et al. (J2001 1. The 
earlier 12 ks RGS observation is too short to include useful high 
resolution spectral information. As the warm absorber could be 




O - 



o 



8 9 10 11 12 13 14 

Observed wavelength (A) 

Fig. 8. Detail of the RGS spectrum between 8 and 14 A. The 
solid line through the data points is our best fit model as de- 
scribed in the text. RGS 1 and 2 are added for clarity in this 
and the next two spectra. 



variable on timescales much shorter than the time difference 
between both observations, we did not use the earlier RGS data 
in our analysis. There are only a few absorption lines that can 
be easily identified. We model the absorber using a combina- 
tion of hot and xabs components. The hot model calculates the 
transmission (both lines and continuum) of a thin slab in col- 
lisional ionization equilibrium. The xabs model does the same 
for a photoionized model. In both models the total column den- 
sity, outflow velocity v and velocity dispersion (Gaussian a) 
a v are free parameters, as well as the temperature T for the 
hot model and the ionization parameter £ = Llnr 2 for the xabs 
model. For more details see Kaastra et al. (2002a 2002b I. To fit 
the spectrum self-consistently we need at least seven absorbers 
(see Table[5}. The best fit spectrum including these absorption 
components is shown in Fig [8] — [K)] Component 1 is the ab- 
sorption due to our own Galaxy discussed in Sect. 14.11 The 
absorber with the largest opacity is neutral and due to the host 
galaxy of IC 4329 A, component 2, and is easily identified from 
its redshifted O I absorption line (see Fig. I10i. Fitting this ab- 
sorber improves the x 2 from 100398 to 16872 for 1603 de- 
grees of freedom (dof) for the continuum parameters quoted in 
Table|2] This absorber is discussed in Sect. 15. II 

The absorber with the highest column density, component 
5, is highly ionized, and including this absorber further reduces 
X 2 to 13559 for 1607 dof. This component is responsible for the 
O VII and Ne IX resonance lines, the N VII Lya and O VIII Lya 
lines; as well as absorption edges due to these ions. Further, a 
less ionized absorber, component 4 with also a lower column 
density is necessary to fit the O V — O VII lines and some iron 
M shell Unresolved Transition Array (UTA) (Behar, Sako & 
Kahn 2001(. Adding component 4 reduces the x 2 to 1 1947 for 
1611 dof. A very highly ionized absorber could also be present 
(component 6), but is only weakly detected, and only improves 
the x 2 by 507 for 4 extra degrees of freedom. This component 
produces (a part of) O VIII and Ne X Lya lines. The fit improves 
significantly, reducing x 2 to 2428 for 1615 dof, by adding a 
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Table 4. Comparison of the Fe Ka emission line as observed in earlier observations. 





ASCA" 


RXTE" 


BeppoSAX 6 


XMM-Newton c 


Chandra d 




1997 Aug. 7-16 


1997 Aug. 7-16 


1998 Jan. 02 


2001 Jan. 31 


2001 Aug. 26 


Energy (keV) 


6.37±0.06 


6.3±0.1 


6.5±0.2 


6.5±0.2 


6.30±0.08 


FWHM (eV) 


920±240 


1150±450 


850±590 


410±500 


440±190 


flux e 


0.7±0.2 / 


l.2±02 f 


1.6±0.7 


<0.9 9 


1.9±0.8 


EW (eV) 


180±50 


240±45 


109±44 


<72 9 


110±43 



Done et al. ( 2000 1, the ASCA and RXTE data were taken simultaneous. 
b Perola et al. (19991 
c based on our own analysis of the pn data. 
d McKernan & Yaqoob 1 2004 1 for the fit with two Gaussian lines. 
e inph m -2 s _1 . 
* Calculated from the given EW and continuum level. 
9 The 3 a upper limit. 




16 17 18 

Observed wavelength (A) 



Fig. 9. Detail of the RGS spectrum between 14 and 20 A. The 
solid line through the data points is our best fit model as de- 
scribed in the text. The O VII absorption at z = is indicated, 
but not fit in the plot. 



22 23 24 

Observed wavelength (A) 

Fig. 10. Detail of the RGS spectrum between 20 and 26 A. The 
solid line through the data points is our best fit model as de- 
scribed in the text. The O VII absorption line at z = and the 
O VII forbidden line are indicated, but not fit in the plot. The 
feature labeled as Unid. is further discussed in the text. 



near neutral absorber (component 3). This fits a part of the UTA 
as well as absorption from O II. However, this component is not 
significant if we assume that the gas in the host galaxy is lowly 
ionized. To distinguish between a neutral or a cold absorber 
the detection of O II is crucial. All these absorbers are further 
discussed in Sect. 15.21 

Finally, we detect a weak O VII absorption component (7) 
from z = plasma, which is discussed in Sect. 15. 31 

5. 1. Absorption from the host galaxy: component 2 

This is the absorption component also detected in the optical 
and is responsible for the reddening observed in the optical and 
the UV. The X-ray O I absorption lines are well explained by a 
neutral absorber. This absorption system is the dominant com- 
ponent (2), in the sense of the opacity in the X-ray spectrum. 

To ascertain that this absorber is neutral and not lowly ion- 
ized, we plot the transmission for a neutral and for a lowly ion- 
ized absorber in Fig. The absorption lines from O II around 



22 A and the deep O II absorption line at about 23.3 A are the 
only detectable differences between these transmission mod- 
els. Juett, Schulz & Chakrabarty (2004i determined from a 
study of the interstellar medium with the Chandra HETG that 
the Oil line wavelength is 23.33 A. In the RGS spectrum of 
IC 4329A there is a deep line at 23.37 A ± 0.01 A rest wave- 
length (23.74 A observed wavelength, identified as Unid. in 
Fig. UOt for which there is no straightforward identification. 
The 0.04 A difference is larger than the 8 mA uncertainty in 
the absolute calibration of the RGS, and is inconsistent with 
the fact that we detect O I at the correct wavelength. 

If this line at 23.37 A is identified as an O II line, then the 
absorber of the host galaxy should be redshifted by 500 km s _1 
and not neutral but very lowly ionized. Inflow is not expected 
if the gas is part of the host galaxy. Therefore we conclude 
that component 2 does not contain a significant amount of O II 
and in the further analysis we will assume that this absorber is 
neutral. 
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Table 5. Absorption components in the IC 4329A spectrum as determined from the RGS spectra. The second order was also 
included in the fit, so as to have a more accurate measurement of the velocity width as well as the outflow velocity. The outflow 
velocity was frozen to km s _1 for the component 2 and 3 as some of the stronger lines in the spectrum have uncertain 
wavelength, or are detected only in the form of blends. 



comp. 


1 


2 


3 


4 


5 


6 


7 


origin 


Gal. 


IC 4329A 


IC 4329A 


IC 4329A 


IC 4329A 


IC 4329A 


Gal/ISM 


model 


hot 


hot 


xabs 


xabs 


xabs 


xabs 


see text 


iVg 


0.455 e 


1.732 ± 0.009 


1.32 ±0.03 


0.32 ± 0.03 


6.6 ± 0.4 


2.0 ± 0.5 




log? 6 






-1.37 ±0.06 


0.56 ± 0.10 


1.92 ±0.03 


2.70 ± 0.06 




rj-iC 


0.5 e 


0.5 e 
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Fig. 11. The difference in transmission for the near neutral ab- 
sorber (T = 1.2 e V, thin line) and a neutral absorber (T = 0.5 
eV, dotted line) for the wavelength range between 21 and 25 A. 
The absorption lines in both models are due to O I. The absorp- 
tion lines only detected for the near neutral absorber are due to 
O II; note the deep O II line at 23.3 A. 



5.2. Warm absorbers 

To explain most of the strongest lines we need an ionized 
absorber. The dominant component 5 is modeled by a xabs 
component and has an ionization parameter log £ = 1 .92, see 
Table|5] £ is measured in units of 10 -9 W m. Due to the low 
signal-to-noise ratio we froze the elemental abundances to so- 
lar while fitting the spectrum. This absorber fits the O VII and 
Ne IX resonance lines as well as the N VII Lya and O VIII Lya 
and (a part of) some weaker absorption lines such as Ne X and 
highly ionized iron. There is thus a warm absorber in IC 4329A 
similar to those detected in other narrow-line Seyfert 1 galax- 
ies. The total hydrogen column density detected for this ab- 
sorber ranges from 1.32 x 10 25 m~ 2 for the lowest ionized 



Table 6. The model column densities of the most prominent 
ions in the spectrum of IC 4329A. For each ion we indicated the 
absorption component that is dominant in producing absorption 
from this ion. 



Ion 


log N H 

(m 2 ) 


comp 


Ion 


log N H 
(m 2 ) 


comp 


N VII 


22.1 


5 


Nex 


21.6 


5 


Oi 


22.2 


2 


Feiv 


20.0 


3 


On 


20.6 


3 


Fe VI 


20.1 


3 


Om 


21.8 


3 


Fe xvin 


20.0 


5 


Oiv 


21.6 


3 


FeXIX 


20.0 


6 


Ov 


20.9 


3,4 


FeXX 


20.2 


6 


Ovi 


21.0 


4 


Fe XXI 


20.2 


6 


Ovn 


21.7 


5 


Fe XXII 


20.4 


6 


Ovin 


22.4 


5 


Fe xxiii 


20.3 


6 


NelX 


21.4 


5 









absorber component 3 to 6.6 x 10 25 m~ 2 for component 5. 
These total hydrogen column densities are similar to the hydro- 
gen column density detected in other Seyfert 1 galaxies such as 
NGC 5548 (Kaastra et al. l2002ai Steenbrugge et al. l2003l 

There are some remaining absorption features in the spec- 
trum that can be fit by a relatively cool absorber (component 3, 
l°g £ = — 1 -37), which in particular produces lines and an edge 
of O III as well as some continuum depression due to the UTA 
of iron. We modeled this absorber with a xabs model. To test 
that this lowly ionized absorber is not part of the host galaxy 
we tried fitting this absorber with a hot component. This wors- 
ens x 2 by 21 for 2377 dof, letting the other absorbers free in 
the fit. However, as the difference in the model is mainly the 
continuum depression due to the UTA, no strong conclusions 
can be drawn as to whether this absorber is photoionized and 
is part of the warm absorber or is collisionally ionized and be- 
longs to the host galaxy. Fig. ^] details the models of these 
seven absorption components and the net absorption model. 
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Fig. 12. The transmission of the neutral absorber (dotted line), 
the neutral and lowly ionized absorber (thin line) and the trans- 
mission also including the medium and highly ionized ab- 
sorbers (thick line). 



In most warm absorbers known up to now, a range of ion- 
ization parameters is detected. Adding component 4 (log £ = 
0.56) improves \ 2 by 52 for 2373 dof. Finally, we detect a 
highly ionized absorber with log £ = 2.70, but with a small 
total hydrogen column density. 

No strong narrow emission lines are detected in the soft part 
of the spectrum. We found an upper limit of 1 x 10 50 photons 
s _1 or an EW < 0.6 A for the O VII forbidden line; and 1.4 x 
10 49 photons s _1 or an EW < 0.03 A for the Ne IX forbidden 
line. There are 1 a significant emission features at 8.54 A and 
11.65 A, however, as there is no straightforward identification 
for these features these are probably due to hot pixels or noise. 

5.3. Absorption at z = 0: component 7 

The last absorbing system detected is due to hot gas at z = 0, 
previously detected in the spectra of other AGN (Rasmussen et 
al. l2U0H N icastro et al. EDoU Steenbrugge et al. l2U0H Behar 
et al. 2003). In the IC 4329A spectrum there are two absorption 
lines due to this local gas, i.e. with zero redshift. We identified 
both as due to O VII, the strongest line being the resonance line 
at 21.6 A (see Fig. 1101. For this line we measure an equivalent 
width of 0.03 ± 0.02 A. The absorption line at 18.62 A (see 
Fig. [9) is expected to be weaker and has an EW of 0.01 ± 
0.01 A. Interestingly, no O VIII Lya line is detected. The Ovi 
line is too blended with a redshifted O I absorption line to be 
detectable. The Ov line at 22.374 A is also not detected at z 
= 0. Assuming that the absorbing gas is in collisional equilib- 
rium, we determine provisionally a temperature range between 
0.03 and 0.2 keV. However, from detecting only two O VII lines 
we cannot conclude that the gas is in collisional equilibrium. 
Assuming that the gas is photoionized, which is less likely as 
the source of radiation is not known, we find an ionization pa- 
rameter of log £ ~ 0.8. The presence of these absorption lines 
in the spectrum indicates that using AGN spectra it is possible 



VII VII VI V IV III I II 




21.5 22 22.5 23 23.5 24 24.5 



Observed Wavelength (A) 

Fig. 13. Part of the spectrum detailing the transmission for the 
seven different absorption components. The thick line indicates 
the best fit spectral model, i. e. with all seven absorption models 
applied. The dotted line at the top with only an O VII absorp- 
tion line is component 7, the other dotted line is component 1, 
or Galactic absorption. The solid line on the top represent com- 
ponent 6, the other components are labeled in the plot. In ro- 
man numerals the ionization state of oxygen is indicated, from 
neutral oxygen (I) to O VII (VII). 

to obtain a measure for the amount of highly ionized local gas, 
as well as an indication of its temperature or ionization balance. 

5.4. Broad emission lines 

Compared to the best fit with our absorption model, there ap- 
pears to be an excess at the wavelengths around N VII Lya and 
O VIII Lya in both the RGS and pn spectra. These excesses can 
be fit by either broad Gaussians or by diskline profiles. For the 
diskline profile we used one set of parameters for the inner and 
outer edge of the accretion disk, the emissivity slope q (as de- 
fined in Sect. 14. 2> and the inclination angle i. The result of the 
model with and without these diskline profiles can be seen in 
Fig- El Leaving the wavelength of the lines frozen to the rest 
wavelength in IC4329A, the diskline model gives a substantial 
better fit, lowering \ 2 by 368 for 2367 dof. However even in the 
fit with the disklines, there is still an excess, most pronounced 
at shorter wavelengths. For the Nvn Lya emission line this 
excess produces a narrow emission line like feature at 25.05 A. 

Table compares the model with Gaussian and disklines 
profiles. Adding the disklines reduces the modified black body 
temperature, as the excess is partly fitted by this continuum 
component. However, the newly derived temperature for the 
modified black body is still quite high, compared to other 
Seyfert 1 galaxies. We did not include a C VI Lya emission line 
as the absorption at the longer wavelength is too strong. Note 
that the inclination determined for the disklines is rather high, 
similar to the inclination of the host galaxy. This is in contrast 
to the inclination determined from the hard X-ray Fe Ka and 
Fe XXVI Lya emission lines, which however are also very well 
fit with narrow Gaussians. 
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Fig. 14. The best fit model with the disklines fit (thick line) 
versus the model without these excesses fit (thin line). 

Table 7. Best fit results for the emission lines observed in the 
soft X-ray part of the spectrum. The lines were fitted using 
Gaussians (columns 2 and 3) and disklines (column 4 and 5), 
in both cases the energy of the line was frozen to its rest-energy 
at the redshift of IC 4329A. The disklines parameters were fit 
for both disklines simultaneously. 





Gaussian 


Gaussian 


diskline 


diskline 


norm" 


25 ± 1 


1.4 ±0.2 


20.4 ± 0.5 


4.2 ± 0.3 


\ b 


24.78 


18.97 


24.78 


18.97 


FWHM C 


4.5 ±0.1 


4.1 ±0.3 






EW (A) 


5.8 ±0.2 


0.57 ± 0.08 


3.8 ±0.1 


0.26 ± 0.02 


Rfn 








1.3 


Rout 








400 


q 






1.9 ± 0.4 


i (deg) 






80. 


9 ±0.4 



a in 10 51 photons s . 
6 in A and in the restframe. 
c inA. 
d in GM/c 2 . 



6. Discussion 

We detect no variability in the softness ratio, although the lu- 
minosity varied by about 17 %. Perola et al. ( 1999) observed in 
their broad band BeppoSAX spectrum significant variations in 
the lightcurve between 0. 1 — 100 keV. However, only marginal 
evidence for variations in the hardness ratios were detected. 
Singh, Rao & Vahia ( fTWTl found from EXOSAT data that 
IC 4329A was continuously variable, with one 12% change in 
luminosity over a 20 ks period, corresponding to a characteris- 
tic timescale r = 170 ks. This characteristic timescale is very 
similar to those measured from our lightcurve. They did not 
record variability on timescales shorter than ~ 10 ks. 

In the IC 4329A spectrum there are seven different absorb- 
ing systems, and disentangling them is rather complicated. First 
there is absorption due to cold gas in our own Galaxy (compo- 



nent 1) and warm gas in its surrounding (component 7). On the 
other hand there is strong absorption due to neutral gas in the 
host galaxy, as determined from the reddening measured in the 
optical, UV bands and our X-ray spectrum (component 2). This 
is consistent with the dusty luke-warm absorbers measured in 
NGC 3227 and Ark 564 (Crenshaw & Kraemer l200Tl Further, 
there is the warm absorber, probably originating close to the 
central source. This absorber is photoionized and represented 
by components 3 — 6. 

This warm absorber has a range in ionization parame- 
ters and total hydrogen column density similar to NGC 5548 
(Kaastra et al. I2002al Steenbrugge et al. l2003t . NGC 3783 
(Blus tin et a l. l200l Kaspi et al. l2002l and NGC 7469 (Blustin 
et al. 2003 1. As such, the warm absorber seems very similar 
to other sources observed, with an increase in column den- 
sity with ionization parameter. The outflow velocity measured 
is similar to the lowest outflow velocity (component 5) of 
NGC 5548 (Crenshaw et al.[2003 1, and is smaller than the out- 
flow velocity observed in NGC 3783 (Blustin et al. l2002> . From 
Table[5]there is a slight trend for increasing velocity dispersion 
with increasing ionization parameter from < 50 km s _1 for 
lowly ionized component 3 to 140 ± 75 km s _1 for the highest 
ionized component 6. No trend of velocity dispersion versus 
ionization parameter was observed for the high signal to noise 
ratio LETGS observations of NGC 5548 (Steenbrugge et al. 
2003 ) or other observed Seyfert 1 galaxies. The opposite trend 
was observed in the Seyfert 2 galaxy NGC 1068 by Brinkman 
et al. ( 2002 1 where the lowest ionization parameters had the 
largest velocity dispersion. 

There is a broad excess measured at some wavelength 
ranges, most notably near the O VIII Lya (18.969 A) and N VII 
Lya (24.781 A) emission lines. These are best fit with disklines 
(solid line in Fig.0 and the thick line in Fig. ll4l >. although even 
then there remains excess emission for both wavelength bands. 
The determination of the inclination angle of the accretion disk 
is an important discriminator between the different models ex- 
plaining the narrow absorption lines. Assuming the excesses 
have a diskline profile, we can determine the inclination angle 
(Laor [T99"Tl . namely 80.9° ± 0.4, consistent with the 80.79° in- 
clination angle measured for the host galaxy (Keel 1980 1. The 
emissivity slope is 1.9, similar to previously measured emissiv- 
ity slopes. A Gaussian fit to these excesses is notably poorer, 
but from the UV band we know that broadened emission lines 
can be quite complex requiring several Gaussians to fit one line 
(Crenshaw et al. 2003 1. Further is the FW0I measured for the 
soft X-ray lines similar to the 13,000 km s _1 as measured in 
the optical. It is thus unlikely that these excesses have diskline 
profiles. 

In the EPIC spectra we detect narrow Fe Ka and Fe XXVI 
Lya emission lines. The intensities of these lines are consis- 
tent with the intensities measured with the earlier Chandra ob- 
servation. The intensity of the Fe Ka emission line is consis- 
tent between all observations since August 1997. This constant 
intensity and the narrowness of the line leads us to conclude 
that the line is probably formed at a large distance from the 
black hole. Assuming that the Fe Ka and Fe XXVI emission 
lines have diskline profiles we find a very low inclination angle 
of 20° and a low emissivity slope of about 1, consistent with 
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McKernan & Yaqoob (2004). The resulting profile is very sim- 
ilar to a Gaussian (dotted line in Fig.0. As a result we cannot 
distinguish from our spectra whether these lines have a diskline 
profile or not. As both determinations of the inclination angle 
discussed are inconsistent, we conclude that it is unlikely that 
disklines are detected in IC 4329A with current instrumenta- 
tion. 

7. Summary 

IC 4329A has a heavily absorbed X-ray spectrum. In the com- 
plex high spectral resolution data studied, we detect seven dis- 
tinct absorbers. From IC 4329A we detect absorption from the 
host galaxy, which is neutral; and absorption from the warm 
absorber closer to the nucleus. The warm absorber is similar 
to the warm absorber detected in NGC 5548, and has four dif- 
ferent ionization components. The ionization parameters of the 
warm absorber components span at least two orders of magni- 
tude. Similar to other Seyfert 1 galaxies, most of the gas in the 
warm absorber is highly ionized. We conclude that IC 4329A 
does not have a luke-warm absorber as suggested from a com- 
parative study in the UV band. The lowest ionized component 
modeled here as absorption from the warm absorber could al- 
ternatively be also absorption from the host galaxy. If this is the 
case, then the gas in the host galaxy is lowly ionized instead of 
neutral. Two of the seven detected absorption components are 
not related to IC 4329A: the Galactic absorption and a moder- 
ately ionized absorber at redshift zero. 

In our best fit to the data we need two broadened lines to 
fit the emission of the O VIII and N VII Lya lines. The fit with 
disklines is statistically better than a fit with two Gaussian lines. 
However, from broad emission line studies in the optical and 
UV band, we know that the line profile of these broad lines is 
complex, and poorly reproduced by a Gaussian line. Therefore, 
we conclude that the broadened lines are similar to the broad 
emission lines detected in the optical and UV band, and are not 
related to the disklines detected in some AGN. 

We detect a narrow Fe Ka line, and conclude from the lack 
of variability since the 1997 ASCA observation, that this line is 
probably formed at a distance of several pc from the nucleus. If 
we fit the Fe Ka line with a diskline, we find parameters similar 
to those obtained from the HEG data (McKernan & Yaqoob 
2004 1. However, this profile can not be distinguished from a 
Gaussian line with current instrumentation, and we thus prefer 
the fit with a Gaussian line. 

From the lightcurve we find that the flux varied by 17%, and 
that the characteristic timescale for variability was between 120 
ks and 190 ks. This is similar to previous measured characteris- 
tic timescales in this source. Consistent with earlier results, we 
did not detect variability in the softness ratio for this source, 
indicating that the soft and hard X-ray emission varied simul- 
taneously. 
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